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Abstract 

We investigate the phenomenology of the gravitino dark matter scenario with a 
stau as the next-to-lightest supersymmetric particle at the LHC. For a wide range 
of gravitino masses the lighter stau is stable on the scale of a detector and gives 
rise to a prominent signature as a "slow muon". The direct stau production via the 
Drell-Yan process is always present and independent of the mass spectrum of the 
other superparticles, thus providing a lower bound for the discovery potential of this 
scenario. Performing a careful analysis with particular emphasis on the criteria for 
observing stau pairs and for distinguishing them from the background, we find that 
the 14TeV run of the LHC has a promising potential for finding long-lived staus 
from Drell-Yan production up to very large stau masses. 



1 Introduction 



Supersymmetry (SUSY) with conserved R-parity and a gravitino as the lightest superpar- 
ticle (LSP) is a viable alternative to the most widely studied scenario with a neutralino 
LSP. A stable gravitino is a perfectly good dark matter candidate [U [2] and may even be 
regarded as favored, since it alleviates the cosmological gravitino problem, allowing for a 
higher reheating temperature after inflation OS]- As the superpartner of the graviton, 
the gravitino takes part only in the gravitational interaction. Therefore, the next-to-LSP 
(NLSP) is typically quite long-livedQ For a charged NLSP, this leads to a spectacular 
signature at colliders, charged tracks leaving the detector and no missing transverse en- 
ergy. It could even be possible to capture NLSPs and to study their decays in detail, 
thus measuring the strength of their coupling to the LSP and the LSP's spin [7J. In this 
way, observations of the NLSP could lead to an indirect determination of the nature of 
the LSP 

In this work, we consider a charged slepton NLSP. In the following, we refer to the 
lightest charged slepton as the stau t\ and allow for mixing between tr and Tl, the 
superpartners of the right- and left-handed tau, respectively. Of course, the results are 
also valid for a selectron or smuon NLSP. For a wide range of gravitino and stau masses, 
the stau NLSP lifetime 

is larger than about 10~ 7 s. Then the stau is metastable, i.e., it usually leaves an LHC 
detector before decaying. 

In this scenario catalyzed big bang nucleosynthesis [H] leads to an upper bound of 
roughly 10 4 s on the stau lifetime. While this bound can in principle be satisfied by 
lowering the gravitino mass sufficiently, a short lifetime is also obtained for gravitino 
masses in the GeV range and a relatively heavy stau [HI [10] due to the dependence of 
on m~ . 

Tl 

Previous studies of the LHC phenomenology of metastable staus have concentrated 
on the production via decays of heavier superparticles, assuming specific scenarios for 
SUSY breaking QU E2 [El M, CGS M El M EH EU ED EZJ E3E2E21E6]. See also [27] 
for a comprehensive review of the topic. Here we do not restrict ourselves to a specific 
SUSY-breaking scenario and we focus on the direct production of staus via the neutral 
current Drell-Yan (DY) process, which possesses interesting properties: 

• The DY contribution is independent of all MSSM parameters except and the 
stau mixing angle 9^, enabling a model-independent analysis. 

• The DY process is always present. Together with the previous point, conservatively 
this leads to an assured discovery potential and strict exclusion limits in a class of 
scenarios characterized only by and 9^. 



1 The same is true in scenarios with an axino LSP, whose interactions are strongly suppressed by the 
large Peccei-Quinn scale, see, e.g., [5j. The NLSP can also be long-lived if its mass is very close to that 
of a neutralino LSP [Bj. We do not study these alternatives in detail but expect the same results as in 
the case considered. 
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Thus, it is natural to ask for the required luminosity which provides this robust exclusion 
limit (and the discovery reach) in the 7TeV and 14TeV LHC run. Although the DY 
production of staus has been included in some studies )28[ [29] 126] . the focus of these 
works has been different and — to our knowledge — this question has only been addressed 
in a brief remark in the review |30j . where the aim is only a rough estimate with fairly 
conservative assumptions. Here, we perform a careful analysis, in particular examining 
the dependence on the imposed cuts and using the proper statistics for small event 
numbers. We also take into account the latest information from the LHC experiments on 
discriminating heavy stable charged particles from muons. We find that the opposite-sign 
stau pair from DY production allows for a clean signal region up to very large integrated 
luminosities. Thus, in spite of its small cross section the DY process is able to provide an 
interesting discovery and exclusion potential, even for relatively large stau masses and 
significantly better than estimated in |30| . 

The paper is organized as follows. In section [2] we examine the DY process, discussing 
the criteria for observing stau pairs and for distinguishing them from the background. 
This allows us to derive the LHC's discovery reach and exclusion potential in section[3j In 
section H] we briefly compare DY production with the production of staus from the decay 
of heavier superparticles. The purpose of this consideration is to estimate for which SUSY 
spectra the exclusion limit from direct DY production is tight and for which spectra it 
tends to be overly loose. We will see that the direct DY production can be dominant for 
large mass gaps between the stau NLSP and the colored superparticles. Unless noted 
otherwise, we discuss the 14TeV LHC run in what follows. 

2 Staus from the Drell-Yan process 

Let us first discuss in detail the DY production of staus, the expected signal in a detec- 
tor at the LHC and the background suppression. All events have been generated with 
MadGraph/MadEvent 4 [31] and its Pythia [32] interface. For the event generation 
in MADEVENT and PYTHIA we enabled the MLM matching scheme [33]. The CTEQ6L1 
PDF set |34| has been used. 

2.1 Background 

Staus as heavy metastable charged particles usually leave the detector. This leads to a 
signal in the tracker and muon chambers. Muons are the only background. Therefore, 
we first study this background in order to devise suitable cuts, which we can then take 
into account in the calculation of the signal in the next subsection. 

The di-muon rate for the moderate and high p^-iange is much smaller than the single 
muon rate and the possible sources are considerably fewer in the case of di-muons. For 
a px-cut smaller and around 15 GeV, b- and c-decays are the dominant sources of di- 
mouns. Above 15 GeV the DY production begins to dominate |35| . In addition to a high 
pr-cut, the b- and c-contributions can be further reduced by isolation cuts. This feature 
relies on the fact that b and c quarks are always produced close to jets while muons 
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from heavy mother particles (like Z or t) tend to be well-separated from the other decay 
products — they are isolated [36J. 

As preselection cuts on the data, we require two opposite-sign muon-like particles 
each satisfying 

• pr-cut: px > 50GeV 

• Barrel cut on the pseudorapidity: \rj\ < 2.5 

• Isolation cut: AlZ^jet = \J Arf^ - et + - et > 0.5 for jets with pp > 50GeV. 

With the px _ an( i the isolation cut, a sufficient rejection of the b- and c-contributions 
should easily be obtained. This is why we refrain from running a detector simulation for 
this issue. All the same we made no effort to specify the isolation algorithm. Instead 
we apply the respective cuts directly on the remaining background (and signal) at the 
level of the PYTHIA Les Houches Event output. The pp-cut also rejects muon pairs from 
on-shell Z decay. Therefore, an additional cut on the invariant mass of the muon pair, 
which was used in [28], would not have a significant effect. 

As the remaining background, we consider di-muons from the DY process and from 
ti-production. We include contributions up to order (two jets) in the case of the DY 
process and (three jets) in the case of it-production and generate 2 x 10 5 unweighted 
events for the analysis. The total normalization of the cross section was fixed from the 
leading order DY di-muon production (without jets) from MadEvent multiplied by a 
constant i^-factor, conservatively chosen to be 1.4, to account for next-to- leading order 
(NLO) corrections. The total di-muon cross section after applying the above preselection 
cuts is then o"b — 25 pb. 

Velocity measurement 

The crucial tool for distinguishing muons from staus is the velocity measurement. A 
significant fraction of staus with a mass of several hundred GeV will have a velocity well 
below the speed of light, whereas the background muons are always ultra-relativistic. 
Unfortunately, measuring the velocity is much more involved than measuring, for exam- 
ple, the momentum. Therefore, the experimental uncertainty is much larger, and a cut 
on the velocity will not reject all background muons. 

There are two distinct ways of measuring the velocity of muon-like particles, the ion- 
ization energy loss (dE/dx) inside the tracker and the time-of- flight (ToF) measurement, 
which measures the time between the bunch crossing and the passing through the muon 
chambers. The former measurement only provides information up to f3 < 0.9 while the 
latter is mainly limited from below — at the design luminosity of the LHC, particles with 
velocities less than about 0.6 cannot be assigned to the correct bunch crossing anymore. 

Although the precision of each velocity measurement is not overwhelming, a combi- 
nation of both measurements provides a highly efficient background rejection. This is 
due to the fact that for stau signal events these two measurements are clearly correlated 
while for background events no correlation is present [37] , According to [38] a cut using 
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Figure f: Left: Pseudorapidity distribution of the considered di-muon background after applying the 
preselection cuts, where if 1 denotes the pseudorapidity of [i ± . Right: Total di-muon cross section after 
applying the preselection cuts but as a function of the variable px-cut discussed at the end of section [3] 
(again on both muons). Both plots are valid for the 14TeV LHC. 



both measurements, fidE/dxi PtoF < 0.8, leads to a background rejection factor of about 
lCT 7 for single stau candidates. 

Thus, if the probability of a mis-identification of two muons within one event is not 
correlated, a background rejection factor of rp = 10 -14 could be achievable]! However, 
for rp = 10 -14 , the relevant background appears only at very high luminosities, when 
pile- up from different bunch crossings becomes relevant. This fact might lower the back- 
ground rejection with respect to the naive expectation of rp = 10~ 14 . To our knowledge, 
currently no quantitative study exists about this issue. In any case, we shall see later 
that a sufficient background rejection (which enables a three-event exclusion over the 
whole considered region for m^) can already be achieved with rg ~ 10~ 10 . Thus, using 
this value is both sufficient and conservative. To show the effect of a looser background 
rejection we will also display the results under the pessimistic assumption that we do not 
gain anything by requiring two stau candidates, hence applying rp = 10 -7 . 

Following these considerations, we will apply the cuts 0.6 < (3 < 0.8 on the signal 
events in order to ensure both a working ToF measurement and a sufficient background 
rejection. 

2 A possible source of such a correlation would be the presence of a highly correlated ?7 + -?7 _ -distribution 
together with a strong dependence of the velocity resolution function on r\. However, figure Q] shows that 
already the former is not the case. 
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Figure 2: Left: Ratio between the electromagnetic contribution and the complete electroweak cross 
section © for = as a function of y/I for up-type quarks (red dashed line) and down-type quarks 
(blue dotted line). In each case the ratio is normalized to be one at s — > oo. Right: Cross section for 
direct di-stau production pp — > Z, 7 — > t^t^ as a function of the stau mixing angle. The curves are 
obtained from a MadEvent simulation with = 500 GeV (solid green, normalized to be one at 9f = 0) 
as well as directly from Q for up-type quarks (red dashed line) and down-type quarks (blue dotted line) 
for a center-of-mass energy y/I = 1000 GeV. The two latter curves are absolutely normalized to be equal 
to the MadEvent prediction at their intersection point. 



2.2 Signal 

The direct stau production via DY only involves two parameters of the more than 100 
free MSSM parameters, the stau mass and the mixing angle Of. 

Dependence on the mixing angle 

We define the stau mixing matrix via 
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and < rrif 2 . The dependence of the di-stau cross section on the mixing angle Of can be 
discussed by considering the tree- level parton-level cross section for qq — > Z, 7 — > t^tJ , 
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Here q denotes the flavor of the annihilating quark, while Q q and T q are its electric 
charge and the third component of its weak isospin, respectively. Besides, c\v = cos 9w, 
% = sin#W) an d s, t, u are the Mandelstam variables. (Taking into account the width 
Fz <C Mz is not vital for the following argumentation.) 

A change in 8^ has an impact on g^, only. Thus, it alters the ratio between the 
three terms in square brackets in ([3]), which are the electromagnetic, the weak and the 
interference terms, respectively. This change is almost independent of the kinematics. 
The terms in square brackets contain s as the only kinematic variable, and even this 
dependence becomes negligible when exceeding a few times Mz- The left panel of figure [2] 
shows the ratio between the first term and all terms in ([3]) as a function of s (arbitrarily 
normalized). From a few hundred GeV on the ratio is almost constant. Thus, in this 
region a change in 0? will only shift the overall cross section without any impact on the 
kinematic distributions. In fact, the dependence of the cross section on 8^ shown in the 
right panel of figure [2] is applicable to all stau masses considered in the following. 

For the plots in figures [3] and H] we consider the case 8^ = 0. However, the minimum 
of a ew (9r) is at ^S 1111 ^ and it is about 7% lower than the value at <7 ew (0). Therefore, 
when estimating the discovery potential and exclusion limits we conservatively choose 
6~ = ^E 1111 . The limits for other values of 8? can easily be obtained from the displayed 
curves. Since it turns out that we can achieve a (almost) clean signal region, the required 
luminosity is, to a very good approximation, simply proportional to the inverse of the 
cross section. 

Dependence on the stau mass 

To show the dependence on the stau mass we simulated the DY production in 
21 mass steps from 100 GeV to 1000 GeV, generating 5 x 10 4 events for each one and 
again considering diagrams up to order a 2 (two jets). We obtained the normalization of 
the total cross section from the corresponding leading-order computation (without jets) 
from MadEvent, corrected by a constant K- factor of 1.35 |39] , This value was found 
for > 200 GeV considering NLO QCD corrections. Additionally including SUSY 
QCD contributions at NLO and next-to-leading logarithmic accuracy yields a if-factor 
between roughly 1.29 and 1.36 001 HQ IS]- 

Figure [3] shows the results. The cuts on 77, p^, and AUrjet have a minor impact on 
the data, whereas the velocity cut lowers the signal by about one order of magnitude but 
with a slightly decreasing tendency when going to very high masses due to the increase of 
slower staus. The fraction of events passing the velocity cut is 8% at = 200 GeV and 
about 20% at = 800 GeV (cf. figured! top). This is due to the fact that the parton 
luminosities in a 14 TeV pp collision begin to decrease more drastically for center-of-mass 
(CM) energies above roughly 1 TeV. The top panel of figure|2]also shows that the velocity 
cut at (5 = 0.6 has considerably less impact on the data than the one at /3 = 0.8. 

The partonic process considered in ([3]) favors perpendicular scattering: (ut — rnf) dt oc 
sin 2 8 d£2, where 8 is the angle between the produced staus and the beam axis in the CM 
frame and dfi the corresponding solid angle element. Thus, the very low region is 
suppressed. On the other hand, the faster decrease of the parton luminosity with increas- 
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Figure 3: Cross section for direct di-stau production pp — > Z, 7 — > t^t^ for 6^ = (ti = tr) as 
a function of at the 14TeV LHC. The impact of the preselection cuts (\t]\ < 2.5, pt > 50 GeV, 
AHfj e t > 0.5) and the velocity cut (0.6 < /3 < 0.8) is displayed. 

ing CM energy determines the high-px tail. This leads to a pT-distribution that peaks 
roughly at pi ~ m^, at least for the mass range < 400 GeV. For larger masses the 
behavior of the parton luminosity at high CM energies shifts this peak a bit downwards 
(see figure |H bottom). Figure HJ bottom also shows that the pp-distributions of the two 
staus are clearly correlated. 

3 Discovery potential and exclusion limits 

In the following we are interested in the integrated luminosity Cjfrn^) required to dis- 
cover or exclude the considered scenario characterized by the parameter m^. The ex- 
pectation value for the number of signal events S is given by 



where c T s(?Wr 1 ) is the signal cross section and e is the detector efficiency The expected 
number of background events reads 



where ct-q is the background cross section and rp is the background rejection factor due 
to the velocity discrimination. (Conservatively, we set the detector efficiency for the 
background to one.) 
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Figure 4: Di-stau distributions at the 14 TeV LHC after the preselection cuts for = 200 GeV (left) 
and = 800 GeV (right). Top: Unweighted events as dots in the /3~-/3 + -plane, where /3 is the 
velocity of r . The black dashed square denotes the region selected by the velocity cut. It contains 
1146 (left) and 2824 (right) events. Bottom: Distribution in the p^-p^ -plane of the events passing the 
velocity cut. 
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Figure 5: Integrated luminosity at which a 5<r-discovery (green solid lines) and a 95% C.L. exclusion (red 
dashed lines) of directly produced stau pairs is to be expected. We have chosen the stau mixing angle 
(9™ ln that yields the smallest production cross section. The dependence on the background rejection 
factor for the velocity cut is illustrated by displaying each curve for r = 1CT 10 as well as rp = 10 
The dark and light red-shaded band around the rp — 1CP 7 curve displays the la and 2<r probability 
band, respectively. The results are shown for the 7TeV and 14TeV run of the LHC. 




Since we expect to obtain solutions that involve small event numbers S and B, we 
consider Poisson statistics. A 5er-discovery corresponds to a set of S and B that fulfills 

1 - e ~ fl E ^r = 3x10 " 7 ' ( 8 ) 

n=0 

where 3 x 10~ 7 is the one-sided p-value corresponding to a 5<7-evidence. 
A 95% C.L. exclusion corresponds to S and B satisfying 



P -(B+S) v w 

^ Z^n=0 n! 

B B» 



0.95. (9) 



In contrast to the case of discovery, the additional parameter appears. This is the 
maximum observed event number up to which the 95% C.L. exclusion is demanded to 
hold. For N = B we obtain the central value of the exclusion limit. Repeating the 
analysis for N at the boundaries of the 1<t and 2a intervals around B then yields the la 
and 2a probability bands for the exclusion limit. 

Inserting S and B as functions of Cj and a-Q according to ([75]) and ([7J turns the 
formulae ([8|) and ([9]) into implicit functions determining /^/(m^). 

As signal cross section we use <ts (m^ ) at % = 05 1111 after applying all cuts mentioned 
above. The trigger and reconstruction efficiency for single stau candidates can be conser- 
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vatively estimated to be 80% |43| . thus we choose (again conservatively) e = 0.8 2 . The 
background cross section <tb is the di-muon cross section after the preselection cuts. As 
explained earlier, we consider the values 10~ 10 as well as 10~ 7 for rp. Figure [5] shows 
the luminosity £j(m^) at which a 5<7-discovery can be expected (green solid lines) and 
at which all scenarios with a metastable stau can be excluded at 95%C.L. (red dashed 
lines, with red-shaded regions around the 10 -7 curve indicating the la and 2a probability 
bands) H The borders of the gray-shaded regions denote the luminosities that correspond 
to 1, 10, and 100 events. We see that both discovery and exclusion are expected to occur 
on the basis of very few events. 

Observing no events when three are expected by a hypothesis is sufficient to exclude 
this hypothesis at 95% C.L. Thus, for luminosities that lead to a sufficiently small back- 
ground B <C 1, the exclusion limit corresponds to the three-event line. This is why for 
small the exclusion limits for different rp are degenerate at the three-event line, which 
coincides with the rp = 10 -10 line in figure [5] — the corresponding background is suffi- 
ciently suppressed. As it is not possible to obtain less than zero events, the red-shaded 
regions do not continue beyond the three-event line. 

In figure [5] we also show the results for the 7TeV run of the LHC. The calculation 
has been analogous except that we have used a slightly smaller if-factor of 1.3 for the 
normalization of the stau production cross section. We find that as the integrated lu- 
minosity exceeds an inverse femtobarn, the LHC is close to tightening the LEP bound 
°f m fi ^ GeV, which currently remains the best limit on the direct production of 
long-lived sleptons. If the 7TeV run reaches 10fb _1 , one will be able to exclude stau 
masses up to roughly 170 GeV. 



Optimized p T -cut 

Looking at the pp-distribution of the staus (figure SI bottom) and the p^-cut dependence 
of the di-muon cross section (figure [TJ right) we see that we can improve the search by 
optimizing the p^-cut according to the stau mass hypothesis being considered. In other 
words, we repeat the previous analysis with an additional p^-cut (again on both stau 
candidates) that grows linearly with m^, 

p T > pjp in = 0.4 x m ?1 . (10) 

The stau cross section after applying this additional cut is nearly identical to the lower- 
most curve in figure [3l We emphasize that the result is not sensitive to the choice of the 
factor 0.4. Choosing 0.5 instead gives almost the same result. 

The dependence of the background cross section on this optimized p^-cut is shown in 
the right panel of figure [TJ Apart from lower statistics, in the high-px region we expect a 
larger theoretical uncertainty. Therefore we perform a conservative fit to the data, such 
that the fitted curve lies completely above the simulated curve in the region j>t > 50 GeV. 

After this change, discovery is expected to take place for one to three events in the 
whole considered region, as shown in figure [6j For instance, a scenario with a 300 GeV 

3 Since even the 2a probability band around the 10 -10 curve is almost degenerate with the curve 
itself, we refrain from displaying the bands for rp = 1CP 10 . 
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Figure 6: Same as figure [5] but with an additional cut pt > 0.4 m^. 

stau NLSP is expected to be discovered at about 10 fb" 1 (and even for rg = 10~ 7 at 
about 20fb _1 ) through direct production alone. On the other hand, if it is not chosen 
by nature the same scenario can be excluded at 95%C.L. with roughly 30 fb -1 . In the 
very long term, nearly the whole considered mass range is accessible at the LHC, for 
example, masses up to 600 GeV (exclusion) and about 750 GeV (discovery) for 300 fb -1 . 
Note that the probability bands for the rp = 10 -10 exclusion curve are degenerate with 
the three-event line. In figure [6] only the 2a probability band for the 10 -7 curve is visible. 

Besides, the dependence of the LHC potential on the background rejection factor rp 
can be studied. The large change of rp by three orders of magnitude causes only a small 
change in the luminosity required for a discovery by a factor of about 2. Note that the 
discovery curves are the expected discovery reach. Since for rp = 10~ 10 the sufficient 
event number for discovery is one or two, the statistical fluctuation is of 0(1), too. So a 
discovery can easily take place at half or double the integrated luminosity shown in the 
figuresO Thus, the difference between the rg = 10 -7 and rg = 10 -10 discovery curve 
is of the same order as the expected statistical fluctuation. For this reason we assess 
the discovery curves to be relatively insensitive to the change in rg in the considered 
range. Moreover, the exclusion limits for rp = 10 -7 and 10 -10 become degenerate at the 
three-event line. Thus, the exclusion potential is not sensitive to rp at all. 

This also sheds some light on the impact of the uncertainties in the background cross 
section. The variation of rp over three orders of magnitude reveals how little impact an 
uncertainty of tens or even hundreds of percent in the background cross section would 
have. This justifies the somewhat arbitrarily chosen iT-factor for the muon background 

4 For this reason we refrain from quantizing the required event numbers to integers. 
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in section 12.11 and shows that setting the detector efficiency for the background to one 
was not very conservative. Similarly, it shows that our result does not rely on the exact 
shape of the high pp tail of the distribution in the right panel of figure Q] 

To realize this optimized p^-cut in an experimental analysis, one could choose the 
value of p™ in corresponding to the stau mass that — for example, according to figure [6] — is 
within reach at the luminosity of the analyzed dataset. 

4 Direct production versus other channels 

The cross section for the production of superparticles heavier than the stau can be larger 
than the DY cross section considered so far. Such sparticles promptly decay into the 
NLSP, emitting Standard Model (SM) particles which provide an additional signature 
with potentially higher significance than the detection of the metastable stau. However, 
the SM particle radiation from cascade decays depends strongly on the sparticle mass 
spectrum and has to be distinguished from a higher SM background. Hence, we assume 
that the easiest way to find SUSY in the considered scenario is the direct detection of 
the NLSP, independently of its production channel. Then we can estimate whether staus 
from direct production or those from decays of heavier sparticles will be the dominant 
contribution to the discovery of a gravitino-stau scenario. This enables us to decide 
for which SUSY spectra our previous calculations for the direct DY production yield a 
good approximation for the potential of the LHC and for which spectra they are overly 
conservative. 

We classify three sources of the production of SUSY particles, 

• production of sleptons, including the direct production of staus, 

• production of neutralinos and charginos, and 

• production of colored sparticles. 

Let us look at the leading contributions of each class. For simplicity, in each case we 
consider exemplary production rates of a single sparticle species and do not sum, e.g., over 
the generations of sfermions (which would require an assumption on the relation between 
their masses). It is easy to estimate the production rate in the case of degeneracies by 
multiplying the results by the appropriate factor. For this consideration we computed 
the cross section via MadGraph/MadEvent 5 03] at lowest order and cross-checked 
whether NLO computations \39\ HO) H5] lead to roughly the same conclusions. 

At lowest order in the electroweak coupling, 0(a 2 ), direct production of sleptons is 
only possible via neutral current and charged current DY. According to figure [21 right- 
and left-handed slepton pair production differs by a factor of about 2 in the cross section. 
On the other hand, such an increase of the cross section is compensated already by a 
rather small increase of the slepton mass by a factor of less than about 1.3. Accordingly, 
the contribution from another neutral current DY produced slepton decaying into the 
NLSP will only be noticeable if T\ — tr and for a very small gap between and fnj^. 

Production of Zl v\ via is in principle enhanced relative to Z, 7 — > l£ l£ by a factor 
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of 2 to 3. Since only Zl couples to VT ± , direct production of T\ via is suppressed 
by sin 2 9^ and thus unlikely for t\ — tr. By contrast, production of it Ui via W + can 
be of the same order as Z, 7 — > t^t^ if mj^ ~ < 1.6t71^. (The production via 
W + is slightly enhanced against the one via W~ due to the charge asymmetry in the 
initial state.) However, in any case, the production of another slepton pair V^vi or / + /~ 
decaying into a pair of t\ does not have the potential to exceed the direct DY production 
drastically. 

Production of neutralinos and charginos is accessible at 0(a 2 ) either via neutral 
and charged current DY or via t- and u-channel squark exchange. DY production of 
neutralinos only occurs in the case of a noticeable contribution of higgsinos and even 
then the cross sections are quite small. In contrast, Z, 7 — > W + W~ has quite a large 
cross section. The same is true for — > W°. For m^ ± < 3 m,^ R the DY production 

of W + W~ can exceed Z,j — > r^r^. If the charginos are more higgsino-like the cross 
sections become lower, m^ ± < 2mf R is required for a competing or dominant higgsino 
production. The production via squarks in the t- and n-channel introduces an additional 
dependence on the squark masses. However, the cross sections are roughly the same 
as the corresponding production via DY only if m-g ~ Thus, for larger squark 

masses the t- and li-channel squark contributions become subleading. According to these 
considerations, one should keep in mind the chargino production via DY, which can 
become an important channel, especially in the case of a light wino. 

Let us finally consider the production of colored sparticles. Since the LHC is a 
proton-proton collider, at leading order 0(a 2 ) squark-pair, gluino-pair and squark- gluino 
production each allow a variety of diagrams that include the initial states gg, qg and qq, 
which are the dominant hadronic channels at the LHC for low, middle and high CM 
energies, respectively. Setting all masses equal to 1 TeV, the production of the colored 
sparticles and direct stau-pair production follow roughly the ratio 

111 

: (TttTj : <r~ : a~+ — ~ 1 : — : — : , (11) 

Q9 99 qq r R r R lg 33 1700 » V J 

where we chose q = u (either right- or left-handed), which is the leading contribution. 
Thus, spectra with dominant direct DY production of staus are required to have quite 
a large mass gap between the NLSP and the colored sparticles to compensate this ratio. 
Large mass gaps are typical for gauge |46| HT] or gaugino mediated |48} 09] SUSY break- 
ing, and no-scale models |50j . for example. Furthermore, (jlip shows that gluino-squark 
production is most likely to be the dominant contribution to the production of colored 
sparticles Jl 

Let us close these considerations with two examples. In a scenario with 
?n ?R = 200 GeV, m^ ± = 600 GeV, = 1.4 TeV, m-g = 1.8 TeV 
or, for a larger overall mass scale, 

= 800 GeV , m^ ± = 1.7 TeV , m s = 2.6 TeV , m-g = 3.2 TeV , 



5 As long as the masses of the squarks are not vastly different, this conclusion remains to hold if we 
sum over the squark flavors and vary the overall mass scale within reasonable boundaries. 
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the cross sections for t^t^, W + W~ and ug production are all roughly of the same size. 

However, in the considered case of a rather large mass gap to the colored sparticles, 
the staus from cascade decays have significantly higher velocities than directly produced 
ones (due to the large phase space) |28^ 15 lj . Consequently, staus from cascade decays 
are more likely to be rejected by the velocity cut needed for the discrimination against 
muons. Thus, already in the exemplary equal-production-rate scenarios from above we 
expect direct DY production to be the dominant contribution to detectable staus. 

5 Conclusions 

Metastable charged super symmetric particles lead to prominent signatures in the detec- 
tors at the LHC. We have shown that these signatures enable a very efficient background 
rejection. As a consequence, despite its relatively small cross section direct Drell-Yan 
production of metastable charged sleptons has an interesting potential for discovering or 
excluding their existence at the LHC for a wide range of masses. Above all, it provides a 
robust lower limit on the LHC potential for scenarios with a metastable charged slepton 
that depends only on the slepton mass. For instance, the 7TeV run will improve the 
LEP limit in the near future and could exclude slepton masses up to roughly 160 GeV 
with an integrated luminosity of 10fb _1 . 

Particularly for the heavy mass range, we have proposed an additional cut depending 
on the slepton mass which may further reduce the background. At the 14TeV LHC, 
this would allow to discover a 300 GeV stau, for example, at about 10 fb" 1 via direct 
production. With a very large luminosity of 300 fb _1 , masses up to roughly 600 GeV 
could be excluded, and even heavier staus could be discovered. As mentioned, this 
mass region can be regarded as interesting from a cosmological point of view due to the 
constraints from catalyzed big bang nucleosynthesis. 

In the spirit of a model-independent approach, we have assumed the stau mixing angle 
yielding the minimal production cross section, which is slightly below the cross section for 
a pure tr. The limits on tl are correspondingly tighter due to its larger production cross 
section. Concerning experimental issues, we chose to be conservative as well. The LHC 
potential may improve, for example, with a better control over the distinction between 
charged sleptons and muons or a better detector efficiency than assumed here. 

By considering channels other than direct Drell-Yan production, a larger, albeit more 
model-dependent, discovery potential and tighter exclusion limits can be achieved. We 
have discussed briefly the production processes that are most likely to be the dominant 
ones if the mass gap between the metastable charged slepton and the heavier superpar- 
ticles is not large enough to guarantee dominant Drell-Yan production. At the LHC, for 
many mass spectra this would be wino production or associated squark-gluino produc- 
tion. 
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